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. 1 Introduction 

q ; 

Q\ ■ In this work we discuss the stability of traveling wave (TW) solutions of the following family 

2 [ of convection-reaction-diffusion equations: 

>! TUtt + u t + g{u)u x = [k(u)u x ) x + f(u). (1) 

•rH . 

Here r > 0, f(u) and g(u) are peace-wise continuous functions, k{u) is smooth and positive 
for u > 0. Equations belonging to this family attracted attention of many authors. The 
case t = 0, recognized as convection-reaction-diffusion (CRD) equation, is the subject of 
investigations in several monographs [TJ [31 S] - It is widely used to describe transport 
phenomena in porous media [5] , theory of combustion and detonation [6] , and mathematical 
biology [3|8]. 

Besides the various applications, the CRD equation is valuable as the simplest nonlin- 
ear model of transport phenomena having in some cases nontrivial symmetry and, thus, 
possessing many exact solutions and conservation laws [9j [T0l [TTJ [T2j [131 EH fl5] . 

General concepts leading to the family (JT|) with r ^ can be found in papers [HI [T7J 
ITS! IT9"] . It can be formally introduced if one changes in the balance equation the convenient 
Fick's law 

J{t, x) = -Wq(t, x), 
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stating the thermodynamical flow-force relations, with the Cattaneo's equation 

d J(t, x) T/ . _ , . 
t g t ' + J(t, x) = -Wq(t, x), 

which takes into account the effects of memory. 

Physically meaningful TW solutions to Eq. (JT]), such as periodic, kink-like and, soliton- 
like solutions, compactons, shock fronts, cuspons, and many other are either shown to exist 
or exactly constructed in papers [HI EE E21 1231 HS1 E£3 E25] • 

The aim of this work is to analyze the stability of kink-like and soliton-like TW solutions 
to Eq. ([1}). The structure of the work is following. In section 2 a geometric insight into 
the family of TW solutions is made, and the linearized equations describing the evolution 
of small perturbations of TW solutions, taking the form of spectral problem are derived. In 
section 3 some important properties of the continuous spectrum are stated, and the stability 
of the constant asymptotic solutions is studied. In section 4 some statements concerning the 
discrete spectrum are formulated. In section 5 the results of numerical study of the temporal 
evolution of the kink-like and soliton-like TW solutions, confirming and supplementing the 
results of qualitative studies, are presented. 

2 Geometric insight into the TW solutions, and the 
statement of the problem 

This work is devoted to the analysis of stability the TW solutions having the form 

u(t,x) = U(z), z = x — st, (2) 

where s stands for the velocity of the traveling wave. To begin with, let us make the 
geometric interpretation of some of the TW solutions. Substituting the ansatz (j2J) into the 
source equation, we obtain the following ordinary differential equation: 

AU"(z) = k{U) U' 2 + [s- g(U)) U> + f(U), (3) 

where A = t s 2 — k(U), both symbols (■), and (•)' stand for the derivative. This equation 
can be presented in the form of the dynamical system 

AU'(z) = AW(z), 

A W'{z) = k{U) W 2 + [s - g{U)} W + f{U). 

Now let us formulate a simple statements concerning the dynamical system fll]) and the 
properties of its solutions. 
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Figure 1: Left: the heteroclinic trajectory in the phase plane (U, U'); right: the correspond- 
ing kink- like solution U(z) 




Figure 2: Left: the homoclinic trajectory in the phase plane (U, U')\ right: the corresponding 
soiton-like solution U(z) 



1. The stationary points of the system (jlj) belong to either the horizontal axis, or singular 
line A = 0. 

2. A smooth kink-like solution is represented in the phase plane (U, W) by the hetero- 
clinic trajectory, which does not intersect the singular line A = (see Fig. [1]) 

3. The soliton-like solution is represented by the trajectory bi-asymptotic to a saddle 
point which does not belong to the singular line A = (see Fig. [2]). 

4. The trajectory bi-asymptotic to a saddle point lying in the singular line A = can 
correspond to either compacton, or shock front (fore more detail see [25| 126]). 

Let us note that obtaining the analytical description to solitons, compactons, or shock 
fronts in case of a typical dissipative system like (j3J) is rather difficult. But it is much more 
easy to "capture" the homoclinic trajectory through two bifurcations: the Hopf bifurcation, 
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Figure 3: Birth of the bi-asymptotic trajectory. Left: periodic trajectory created as a result 
of the Hopf bifurcation; middle: the change of driving parameter causes the grows of the 
radius of periodic trajectory, which is finally destroyed as a result of the interaction with the 
nearby saddle point (right) 



followed by the homoclinic bifurcation (see Fig. [3]). The Hopf bifurcation can be predicted 
by means of the local asymptotic analysis [27] [28] . The homoclinic bifurcation is nonlocal, 
and therefore should be captured numerically [29j 122] . 

A simple analysis shows, that the equation (CO) can have the homoclinic or heteroclinic 
TW solutions, if the corresponding dynamical system possess at least two stationary points. 
This, in turn, determines the form of the source term f(u), which in the simplest case is as 
follows: 

f(U) = (U-U 1 )(U-U )^(U). (5) 

We assume that Uq < U±, and ^f(U) does not intersect the horizontal axis within the interval 
(U , U 1 ). 

Let us now concentrate upon the problem of stability of the TW solution ([2]). Since we are 
interested in studying the stability of kink-like and soliton-like solutions, then we assume, 
that 

lim U (z) = m±, (6) 

z— >±oo 

where m± are the constants coinciding with Uq or U\ (in the case of a kink-like solution 
m + 7^ m_, in the case of soliton-like solution m + = m_ = Uq). We assume in addition, 
that 

lim U (k \z) = (7) 

z— >±oo 

for any natural k. 

To study the stability of a TW solution U(z), we use the ansatz 

u(t, x) = U(z) +e exp[\t]v(z), (8) 
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where A is the spectral parameter, and |e| *C 1. It is instructive to pass to new independent 
variables 

t — tj Z — X S t • 

in which the invariant solution (j2j) becomes stationary. In the new variables the equation 
(fH) reads as follows: 



T 







1 2 



dt S dz 



u + 







d 



dt S d z 



u + g{u) 



du d 

dz dz 



du 
dz 



f(u) 



(9) 



(for simplicity, we omit the bars over the independent variables henceforth). Up to 0(e 2 ) 
the function v(z) satisfies the equation 



' d i 

a z 



viz 



+ 



g(U) -2Uk(U) - s(l + 2r A) 



d 

. dz 



+ 



(10) 



+A(1 + r A) + g(U)U - k{U) U - k{U)U 2 - f(U)j v(z) = 0, 
where A = t s 2 — k(U). 

Definition 1. The set of all possible values of A G C for which the variational equation 
( TJ7| ) has nontrivial solutions is called the spectrum of the operator L [z, A] . 

Definition 2. We say that the TW solution U(z) is (linearly) stable, if any possible 
eigenvalue A for which the equation (TjTJ) has nonzero solution satisfies the condition A G 

cu o. 



Remark. It is easily seen, that zero eigenvalue always belongs to the spectrum of the 
operator L \z, -4-, A], for the following statement holds: 

Lemma. If U(z) is a TW solution of the equation (T7]) ; then v u (z) = U(z) is the eigen- 
vector of the operator L [z, 0] . 

Proof. Differentiating ([3]) w.r.t. z, we can rewrite the resulting equation in the form: 



d 



[rs 2 -K{U)] g(U)-2Uk(U)-s 



d 

dz 



+g(U)U - k{U) U - k(U)U 2 - /([/)] v v {z) = 0, 

where v u (z) stands for U(z). The differential operator inside the braces coincides with 
L [z, £, A], when A = 0. □ 

As usually, we distinguish the continuous spectrum a cont C C, and the discrete spectrum 
o~discr C C. Being somewhat informal, we can treat o~ C ont as the subset responsible for the 
stability of the stationary solutions m±, and Odiscr for the stability of the solution U(z) itself. 
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3 Stability of the asymptotic stationary solutions 



Nowwe are going to state the conditions which guarantee that a cont C C~. For this purpose, 
we study the stability of the stationary solution U k = const, assuming that it coincides with 
m + or m_. We assume in addition that the eigenvectors of the operators L [±oo, A] 
belong to the space of tempered distributions S'{1Z) [30]. With this assumption, we can 
solve the spectral problem 



d X 

a z 



v{z) = 0, 



where I stands for plus or minus infinity, by applying to this equation the Fourier transfor- 
mation. Thus, assuming that \\m.\J(z) = U k , we get: 



F { L [4 h ^] «}fcl = 
= L[£, X] F[v](0 = 

-[ts 2 - K(U k )] e - it [9 (U k ) - s (1 + 2 A r)] + A (1 + A r) - / (U k ) \ F[v] (0 = 0. 



(12) 



This equation has nonzero solution F[t>](£) e S'(1Z) if 
-(1 + 2istO±VQ 



X 



2t 



where 



Q = l- Atk? + Arf{U k ) + Air^g{U k 



Since we are going to estimate the real part of A, it is instructive to use the representation 
Q = x + iy, with x, y 6 R. 



Raising both sides of this equation to the second power, and eliminating y, we get the 
bi-quadratic equation 



4 2 
X — X 



Its positive root takes the form 
1 



x 



V2 



l-4r< i + 4r/(f/ fc ) + ^Q 1 



1/2 



where 



Qi = 1 - 4r< 2 + 4r/(l7 fc ) + 16rW(£4)- 



6 



So the biggest real part of A, which we denote by Re\ + , is as follows: 



ReX 



1 

27 



-1 



V2 



1-4t< 2 + 4t/(£4 



1/2 



Let us solve the inequality Re\ + < with respect to £ G 
inequality 

1/2 



It is equivalent to the 



1-47-^+47-/(1/*) + 



which, in turn, can be rewritten as 



/QT < l + 4r< 2 -4r/(f4 



< 



y/2, 



Rasing both sides of this inequality to the second power, we get, after some algebraic ma- 
nipulation, the inequality 



e [rg\U k ) - K (U k )] < -f(U k ), 

which should be fulfilled for any £ e K. So the following statement is true. 
Statementl. The stationary solution U k is stable if 



(13) 



rg\U k ) - K{U k ) < 0, and f{U k ) < 0. 



In order to get a cont in case r = 0, we have to consider the eigenvalue problem 

^ £} v(z) = 

= {<U k )£ + [s- g(U k )] £ + f(U k )) v(z) = Xv(z). 

Applying the Fourier transformation, we get 

\ = ii[g(U k )-s}-K(U k )e + f(U k ), (GR. 

So in the case r = 0, we obtain the following result. 

Statement 2. If for U k = lim U(z) f [U k ] < 0, k = 0, 1 then a cont E C~. 

4 Some remarks concerning the discrete spectrum 

We remind that if the function U(z) is the TW solution of the equation (JT]), then its derivative 
U(z) is the eigenvector of the operator L [z, 0] , corresponding to the eigenvalue A = 0. 
This fact plays very important rule when r = 0, for it is possible to address the question 
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of stability of soliton-like and kink-like solutions by employing the classical Sturm-Liouville 
theory [31] . Let us shortly remind one of its conclusions. 

Theorem (Sturm Oscillation Theorem). Let Xq > A2 > ... be the eigenvalues of the 
spectral problem 



H 



d 

dz 



d 2 

u(z) = I -—7 + Viz) /> u(z) = Xu(z), lim u(z) = lim u(z) = 0, 

1 dz 2 z^A K ' z^B 



where u(z) G L 2 (A, B), (A, B) G R, is finite or infinite interval, V(z) is a bounded func- 
tion. Then the eigenvector u(z, X n ) has exactly n zeroes in (A, B). 

Let us consider the eigenvalue problem for the case r = 0, assuming in addition that 
K = const. Then we can rewrite the variational equation fflUj) in the form 

" d 2 _ d 



+ a(z)— + b(z) 



dz 2 dz 
where A = X/k, 



v(z) = Xv(z), (14) 



aW = i^M, Hz) = m-um 

K K 



The problem (Tl4j) can be presented in the standard Sturm-Liouville form, if we use the 
following transformation: 

v(z) = exp [<f(z)] w(z), where ip(z) = — . (15) 

Using ( fl5i) . we obtain the eigenvalue problem 

~ d 2 



+ $(z) 



^d z 
where 

a 2 (z) a(z) 



w(z) = Xw(z), (16) 



$(2) = b(z) - 



1 



By w u (z) we denote the eigenvector corresponding to the eigenvalue X u = 0. And now, 
if U(z) is the kink with the monotone profile, then v u (z) = U(z), as well as w u (z) = 
U(z) exp [— (f(z)], is a functions, which does not intersect the real axis, see FiglH Then on 
virtue of the Sturm Oscillation Theorem, X v = Xq = 0, and all the remaining eigenvalues are 
negative. 

On the other hand, if U(z) is the soliton-like TW solution, then both v u (z) = U(z), and 
w u (z) = U(z) exp [— ip(zj\ intersect the horizontal axis, see FigEJ So X u = \\ = 0, and there 
is the eigenvalue Aq > belonging to the right half-plane of the complex plane. The latter 
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Figure 4: Graphs of the monotonic kink-like solution (left), and its derivative (right) 



result is rather well-known [32J . Conditions concerning the stability of kink-like solution can 
be formulated as follows. 

Statement 3. A monotonic kink-like solution of the equation (QJ) with r = is stable, 
provided that f(Uk) < for k = 0, 1. 

It is possible to apply the approach based on the Sturm-Liouvliie theory in the case 
when t > 0, and g(u), k(u) are constant functions. Under these conditions, the variational 
equation can be presented in the form 



+ a(z)-£- + b(z) 



dz 2 dz 



v(z) = 0v(z), (17) 



where 



s(l + 2r\)-g ~ f(U) A(l + rA) 

a(z) = , b{z) = , v = 

n n n 

fi = k — t s 2 , The transformation 



v(z) = exp [(f(z) w(z)] , where cp(z) 



g - s(l + 2r A) 



2/x 

leads in this case to the spectral problem 
d 2 



dz* 



w(z) = x(^) w ( z )> hmw(z) = \imw(z) = 0, 

z—i-a z^b 



where 



= g{2s-g)-s 2 + f[U{z)} = \[K-gsT + KT\} 

4 /i 2 A* t^ 

Thus, if U(z) is a monotonic kink-like solution, then Vq(z) = U(z) is the eigenvector 
corresponding to the eigenvalue A = 0. The corresponding function wq(z) = e - ^^ Vq(z) is 
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Figure 5: The graphs of one-humped soliton-like solution U(z) (left), and its derivative U(z) 
(right) 




Figure 6: Graphs of the function x(^) an d the eigenvalues Ag" corresponding to U(z) when 
U(z) is a kink-like TW solution. Left: K — g s r > 0; right: K — g s r < 



the eigenvector of the operator H = + Since the function Wq{z) does not intersect 

the horizontal axis, then, on virtue of the Sturm Oscillation Theorem, it corresponds to 
the eigenvalue xo = 0, and any other eigenvalue of this problem is negative. But x(A) is 
the quadratic function of A, and therefore the source eigenvalue problem (1171) can have an 
extra eigenvalue, corresponding to the function U(z). This extra root will be negative, if the 
constant k — g sr, is positive, and negative otherwise, see Fig. [6j Hence for r > we get an 
extra condition 

K-gsr > 0, (18) 

assuring that Odiscr £ C~ ■ 

If, in turn, U(z) is the soliton-like solution, then the eigenvector w v = exp [— (p(z)] U(z) 
corresponding to the eigenvalue \ v = intersects once the horizontal axis. Hence v — 1, and 
there is an extra eigenvalue xo > Xi = 0, to which corresponds a pair of the eigenvalues A^ 
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Figure 7: Graphs of the function x(A) and the eigenvalues Xq corresponding to U(z) when 
U(z) is a soliton-like TW solution. Left: k — g s r > 0; right: k — g s r < 

of the source eigenvalue problem. As it is seen in Fig. regardless of the sign of k — g s r, 
there exists the positive eigenvalue Xq, hence the soliton-like solution is unstable. 

5 The results of numerical simulations 

In the general case, i.e., when, e.g., the function g(u) is not constant, estimation of adiscr 
is rather more delicate problem. In paper [33] such estimation is performed for the TW 
solution 

W(z) 

u(t, x) = U (z) = — — — — , $f(z) = exp [mi z) + C 2 exp [m 2 z] + C 3 exp [m 3 z], (19) 
satisfying the equation 

Tu tt + u t + [iuu x = ku xx + u(u- mi) (u - m 2 )(u - m 3 ), (20) 
under the following restrictions on the parameters: 

1 + Jl + 4rK(2 =1 m k ) 2 
/i = 3 A, v — — A, s = (m 1 + m 2 + m 3 ) A = ^ 

Unfortunately, cr cont for none of these solutions is contained in C~, as will be shown below. 

Statement 4. // Cf + Cf > 0, and < mi < m 2 < m 3 , then the stationary point 
m + = lim U(z) is unstable. 

z — > +oo 

Proof. The stationary point Uk = m + is stable, if the inequality ( TIB"]) is fulfilled for all 
£ G R. The case C 3 = is very easy to analyze. Indeed, under this assumption m + coincides 
with m 2 , and thus 

f(m 2 ) = v [m 2 — mi) [m 2 — m 3 ) = —A (m 2 — mi) {m 2 — m 3 ) > 0. 
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Hence the conditions of the statement 1 are not fulfilled. 
If C3 7^ 0, then lim = 777.3, and 

z — > +00 

/(^3) — ^ ( m 3 — m l) { m 3 — m 2) — ~ A (^Tl3 — mi) (7TI3 — ?772) < 0. 

Let us address the condition 

t g 2 (7773) — k = t (fims) 2 — k = r(3 A 7773) — k < 0. 
This inequality is equivalent to the following one: 



1 + VI + 4tkM 2 / « w 

A= O ^72 < A/n 2' M = 777i + 777 2 + 777 3 , 

2rM 2 y 9 r 7773 
which can be rewritten as 



Vl + 4 r k M 2 < 2tM 2 J 

y 9 r 7773 

It is evident, that the above inequality cannot be fulfilled if the RHS is negative, so let us 
assume, that 2 r M 2 gT K m2 — 1 > 0. Raising both sides to the second power, we get, after 
some algebraic manipulation, the inequality 



" <-^(M 2 -9ml). (21) 



9 r 777,3 9 m 3 



The inequality (jzTj) cannot be fulfilled since 

M 2 — 9 7773 = (mi + m 2 + m 3 ) 2 - 9 7773 < (3 m 3 ) 2 = 9 m 2 3 = 0. 

□ 

The above statement tells us, that, whether or not cr^ scr belongs to the left half-plane, 
the traveling wave fll9p cannot evolve in a self-similar mode. Example presented in Fig. [H] 
confirms this conclusion. The numerical solution of the Cauchy problem with the Cauchy 
data u(0, x), u t (0, x) being equal, respectively, to \&' \x) /^f(x) and — s (^'(x) /*ff(x)) , per- 
formed under the following values of the parameters r = 1, k — 1, m\ — 0.5, 772,2 = 1-5, 7773 = 
5, C2 = 1, C3 = 3 shows that the initial perturbation evolves for some time in a self-similar 
mode. In the long run the self-similar evolution becomes corrupted. The process of self- 
similarity destruction starts from the far end, and this is in agreement with the observation 
that the condition (Tl3|) is not fulfilled for u = 7773 = const. 

Numerical experiments performed with another values of the parameters, namely: r = 
1, k = 1, 777,1 = 1, 777 2 = 2, 7773 = 3, C 2 = 100, C 3 = 0.01 expose the same tendency, Fig[9J 

Numerical simulations purposed at studying the evolution of soliton-like solutions to the 
equation ([T|) show, that these solutions are rather unstable. In a number of numerical 
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O r x 

5 10 15 20 25 30 

Figure 8: Numerical solution of the system ( 1201) in case when r = k = 1, mi = 0.5, ni2 = 
1.5, m3 = 5, C*2 = 1, C3 = 3, and the invariant kink-like solution f fl9|) is taken as the Cauchy 
data. Successive graphs present the TW, moving from left to right 




20 40 60 80 100 

Figure 9: Numerical solution of the system (1201) in case when r = K — 1, m\ = 1, m 2 = 
2, m 3 = 3, C*2 = 100, C3 = 0.01, and the invariant kink-like solution f fl9|) is taken as the 
Cauchy data. Successive graphs present the TW, moving from left to right 
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r =0, f 1= 8, t 2 =l5, t 3 =23 



u 



1.5 



0.0 



0.5 



1.0 



-25 



-20 



-15 




10 



x 



Figure 10: Numerical solution of the system ( 122]) performed with r = 2, k — 1, 7 = 1, when 
the invariant soliton-like solution obtained by solving the system (123]) is taken as the Cauchy 
data. 

simulations, performed with different functions g(u) and k(u), we had encountered three 
types of instabilities destroying the TW solution. The first type is connected with the 
instability of the constant asymptotic solution. Two other types are manifested in either 
fading or blowing-up of solution in finite time. Since this behavior is rather typical, let us 
illustrate it on the example of the solitary wave solutions of the equation 

TU tt + U t + UU x = ku xx + ju (l-u 2 ) (22) 

considered in paper [20J. The factorized system 



obtained via the substitution (j2J), possesses the homoclinic solutions attained through two 
bifurcations. The first one is the Hopf bifurcation taking place when 7 A > [2T], and the 
parameter s is close to the unity. The second one is non-local and is captured numerically. 

In both possible cases, i.e., when for 7 and A are simultaneously positive or negative, 
variation of velocity s near the unit value leads to the appearance of homoclinic loop corre- 
sponding to the soliton-like solution. 

Numerical integration of the Cauchy problem for (122]) with the soliton-like solution taken 
as the Cauchy data, performed with the following values of the parameters r = 2, n — 1, 7 = 
1, s — 0.84912 (when both A and 7 are positive), reveals that the initial solitary wave evolves 
in a self-similar mode for some time, but finally is destroyed as a result of instability of the 
asymptotic state, Fig. [10J Let us note, that the instability of the stationary asymptotic 



AU'(z) = A W(z), 

A W\z) = [s- U[z}} W + U[z] (1 - U 2 



[*])■ 



(23) 
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1=1,3,5,7,7.8,8.3,9.45 

u 




Figure 11: Numerical solution of the system ( )22|) with r = 0.25, k = 1, 7 = —1 (M = 0.71) 
when the invariant soliton-like solution obtained by solving the system ( 1231) is taken as the 
Cauchy data. 



f =0, f,=6, r 2 =io, h=W.l 

11 

8r 

6 ■ M 
4 - 1 

-15 -10 -5 5 10 * 

Figure 12: Numerical solution of the system (1221) in case when the invariant soliton-like 
solution obtained by solving the system ( )23|) with r — 1, k — 2, 7 = — 1 (M = 1.004) is 
taken as the Cauchy data. 
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solution u = Uq = is confirmed by inspection of the formula (fT3j) . Evolution of the soliton- 
like solution appearing in two other cases depends on the magnitude of the characteristic 
number M = u max /Co, where u max is the maximal amplitude of the initial perturbation, 
Co = a/k/t is the acoustic waves' characteristic velocity. We call M the "Mach number", 
since it plays the analogous role, as the parameter known under this name plays in the theory 
of supersonic flows. Indeed, if M < 1, then the initial solitary wave vanishes to zero, Fig. [TTJ 
Let us note, that we deal in this case with the destruction mechanism, which is completely 
different from the convenient dispersion. The TW evolves for a while in self-similar mode, 
but at some instant its amplitude is subjected to the drastic decrease so that the wave pack 
completely vanishes in finite time. 

When M > 1 then the evolution of the solitary wave ends with the blow-up regime 
appearance, Fig. [T2J We would like to mention in this place, that such a strong dependence 
of solutions of the nonlinear hyperbolic- type equation similar to ([1]) upon the " Mach number" 
was noted for the first time in paper [T8] . 

6 Discussion 

Thus in this work qualitative and numerical investigations of the TW solutions of the system 
([T]) have been performed, with special attention paid to the behavior of the kink-like and 
soliton-like solutions. On the basis of the results obtained, we can state that incorporation 
of the second- derivative with respect to time does not lead to drastic change of the situation 
taking place in the case of the classical convection-reaction-diffusion equation, for which a 
monotonic kink-like solution is stable if both of the constant asymptotic solutions are stable, 
while the soliton-like solution is always unstable. In fact, the situation in the case of the 
equation ([I]) is not so clear, for the Sturm Oscillation Theorem cannot be directly applied 
when g(u) and k(u) are nontrivial functions. Yet in the situation when both k(u) and g(u) 
are constant, and we can formulate the variational problem in the Sturm-Liouville form, 
the results of qualitative analysis show, that an extra inequality should be fulfilled in order 
that the kink-like solution be stable. Under the same conditions, any soliton-like solution 
proves to be unstable, as this is the case when r = 0. Numerical experiments performed 
with different k(u), g(u) and f(u) for which the source equation possesses the solitary- 
wave solutions (taken as the Cauchy data), reveal the instability in the wide range of the 
parameters' values. They evidence that solitons, compactons, and shock waves appearing in 
the class of the generalized convection-reaction-diffusion equations remain unstable in the 
case of positive r. 
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